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Abstract 

The charmonium sum rules for the axial-vector current are considered. The three-loop per- 
turbative corrections, operators up to 8 dimension and ag-corrections to the lowest dimension 
operator are accounted. The contribution of the operators of 6 and 8 dimensions is computed 
in the framework of factorization hypothesis and instanton model. For the value of gluon con- 
densate the following result is obtained: (^G^) = (0.005 + 0.001 - 0.004) GeV^ (for the charm 
quark mass ifi = 1.275 it 0.015 GeV). 
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1 Introduction 



It is well known that the nonperturbative phenomena in QCD are of great impor- 
tance. Indeed, quark and gluon condensates determine the properties of hadrons and 
their interactions to a considerable extent. This was stressed by Shifman, Vainshtein and 
Zakharov in pp, where the technique of QCD sum rules was proposed. 

In particular, authors of jT] pointed out to the significance of the gluon condensate 
in nonperturbative QCD. First of all, gluon condensate has the lowest dimension among 
chirality conserving condensates and that is why it plays dominant role in the sum rules, 
corresponding to the processes without chirality violating. Moreover, the gluon conden- 
sate is directly related to the density of vacuum energy. 

The numerical value of the gluon condensate was found in from the analysis of 
charmonium sum rules. It was 

{-^G%G%) = 0.012 GeV\ 

In this estimation the value Aqcd = lOOMeV was used and a^-corrections of the first 
order were taken into account. 

However, at the present moment Aqcd is known to be sufficiently larger, 
Aqcd ~ 250 MeV, and, furthermore, a^-corrections of the second order to the polarization 
operator are available. In addition, authors of worked at = 0, where the higher 
order terms in the operator expansion series are significant 0. These facts are taken into 
account in the recent paper |3j, where the charmonium sum rules for the vector current 
are considered over again with the purpose to obtain the value of (^G^j^G^j,). (Of course, 
there were many other papers, devoted to the calculation of the gluon condensate, the 
short list of them can be found in [^j, for review see j3]). 

As the result, authors of obtained {^G^^^G"^^) = 0.009 ± 0.007 GeV^ and rfi = 
1.275 ±0.015 GeV for the c-quark mass in MS scheme. Thus one can see that the accu- 
racy of the quark mass value is high, whereas the error in the value of the gluon condensate 
is comparable with its magnitude. 

In the paper j3] the charmonium sum rules for the pseudoscalar current were analysed, 
and for the gluon condensate value the following restriction was obtained: 
{^G%Gp< 0.008 GeV\ 

In the present paper we try to obtain the gluon condensate by considering other inde- 
pendent channel in charmonium, namely, we analyze the axial- vector current. 

It should be mentioned that the values of gluon condensate and c-quark mass are in- 
terdependent in the charmonium sum rule, i.e. the variation of one of them results in the 
changing of another. That is why the determination of independent restriction to each of 
them is a special problem. In the present paper we use the quark mass value, obtained 
in 0, as the input parameter and devote our attention to the gluon condensate. 
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2 Correlator of the axial-vector currents 

We consider the correlator of the charmed quark axial-vector currents: 

I J d'xe^^^T{j,{x)3tm = {q,qu - g.uq'Mq') + q.qMq') ■ (1) 

Here = 07^750, 11^ is the longitudinal part of the polarization operator. In what follows, 
only transverse part of the polarization operator n(g^) is considered. 

n(g^) (1) can be expressed through its imaginary part with the help of dispersion 
relation: 

U(q^) = ^7 ^^^^^^ 
^ 47r2 J sis - g2) ' 

where 

R{s) = 47rlmn(s + iO) , (2) 

m is the pole mass of c-quark. 

The dispersion relation can be saturated by the contributions of physical states. In 
the axial- vector channel the only resonance with the mass = 3510.51 ± 0.12MeV is 
known [H]. We use the simplest model of the spectrum, containing this resonance and 
continuum. As usual, in order to suppress the continuum contribution one should consider 
the derivatives of the polarization operator in Euclidean region = —q^ > 0: 

Usually the quantities M„ are referred to as moments. 
From the phenomenological point of view 



mm? 



(ml + Q2)n 
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and in the ratio of two successive moments the contribution of the lightest state dominates 
at large n: 

where 6c is the continuum contribution. 

The QCD part of the polarization operator consists of the perturbative and nonper- 
turbative terms: 

The first one is determined by its imaginary part (2). R{s) can be expressed as the 
series in the coupling constant a^: 

Ris)= Y: R^''\s,fx')a'ifx'). (4) 

fc=0,l,... 

Hereafter we denote = Q;s(/i^)/7r. We will take first three terms in this series. 

R{s) is the physical quantity and does not depend on the scale /i^, but each term in 
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(4) can depend. First two terms of (4) are known analytically. They do not contain scale 
dependence. One can find them in j7j: 

R^^^=v\ (5) 
rW = 1 lv^({l + v^){2 ln(l - p) Inp + ln(l + p) \np + 2Li2(p) + Li2(p^))- 



-2v(2 ln(l -p) + ln(l + p))) - ^^(21 + 21v + 80v'^ - 16v^ + 3v^ + 3v^) \np + (6) 



1 n'' 

n=l 



In these expressions v is the quark velocity, v = y^l — Aw? / s, p = {1 — v)/{l + v). Lit is 
the polylogarithm function: 

n=l 

The term R?'^ is represented usually as the sum of five gauge invariant parts: 

= clRf + CaCfR^n\ + CpTniRf^ + CpTRf + CpTRf . (7) 

Here Ca = 3, Cp = 4/3, T = 1/2 are group constants and ra^ = 3 is the number of light 
quarks. 

The term R] corresponds to the so-called double-bubble diagram, the diagram with 
two quark loops, in external loop there are massive quarks, whereas in internal one - 
massless quarks. Rf^ has the following form [7]: 

flp) = f-il^^^!-A)fl(') + ^i^ (8) 



4 4s 12 

where 5\ is given by equation (79) in [7]. 

Rp appears from the double-bubble diagram with massive quarks in both loops. In 
the domain s < IQm? only the virtual massive quarks contribute, and has the form 

Rf = 2v'ReF^l--\r.^R^'\ (9) 

The expression for F^ q can be found in appendix B of 0, equation(166). 

For s > 16m^ the contribution of the real quarks appears. It is given by double inte- 
gral, which can not be taken analytically (equation (71) in T). The numerical calculation 
shows that it is small, nevertheless, we will take it into account. 

(2) (2) 

The functions R\ and -R)va appear from diagrams with single quark loop and a num- 
ber of gluon lines, they contain abelian and nonabelian exchanges correspondingly. These 
functions are not known analytically, consequently, one has to use some approximations. 
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Our approximation expressions are based on the first eight moments at = 0, which 
are known analytically [H]. The obtained formulas, being substituted into (3), have to 
reproduce these eight moments with high accuracy. In order to construct such approxi- 
mations we perform the following steps. 

Let us consider the series 

^(2) _ 3 ^ (2) k 
fc=l,2,... 

where z = q^/{Am'^), H^^^* = EH^'^^a'^ and 

k 

n(2) = dlif + CaCfTI-^V + CpTnilif^ + CpTIif + CpTIif similar to (7). The 
coefficients C)^^ d ■••) ^na s known analytically. First of all, we reexpand this series in 
terms of variable uo, 

1 - v^r^ 

Thus we map the complex g^-plane to the unit circle. Then we construct the Fade ap- 
proximation, which usually gives a better accuracy than Tailor series. It has the following 
form: 

ao + aiu + ... + aiUj' 

11 = ; ; . 

^ ^ 1 + biuj + ... + bjUJ^ 

Since we have 8 moments in hand, we can construct Fade approximation with 8 parameters 
Qi and bj. The best results turns out to give the following approximation: 

10.8547 + 9.43221 - 8.76722 cj^ - 1.74256 + 0.853743 cu^ - 0.257734 cu^ 
^ 1 + 0.373686 - 0.439076^^2 ' ^ ^ 

In the similar way we obtain for the abelian part: 



^ 9.25606 - 288.334 cu^ - 394.513 cj^ - 69.8439 cu^ + 19.9673 cj^ - 0.936404 cu^ 

^ 1 - 32.0000 cu- 15.2803 ' ^ ^ 

The last term in (7) is generated by the diagram with two triangle quark loops (singlet 
part). We need the first moments in the expansion 

n?' = T|;j T. cfl^- (12) 

fc=l,2,... 

Such diagrams for different currents were discussed in paper jH]. However, to cancel the 
axial anomaly, in the axial-vector case the current = 07^750 — s^jfj^'j^s was considered. 
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As the result, in the calculated factors Cfj. {k = 1, 7) the logarithms ln(g^/m^) appear 
due to massless cut (see equation (8) in [Hj). For our purposes we need to subtract the 

~ (2) (2) 

contribution of the massless s-quark in Cg f.. The expression for the imaginary part Rg^ 
(the contribution of the massless cuts) can be found in ^Oj- But the corresponding integral 
can not be taken analytically. The purely numerical integration is also impossible because 
of the presence of divergent logarithms ln(g^/m^) (which cancel out after subtracting from 
Cg'l ). The suitable technique was suggested in [S]. Rg'} is represented as the sum 

p{2) _ 9 1 s , p'(2) 
Then the integral over R^J is divided into three parts: 



r — z J r — z J r — z J r — z 

e 1 

In the first integral on the right hand side we expand RgJ near r = and take about 50 

terms to obtain a stable sum of two first integrals in the range e = 0.65. ..0.75 with the 

high accuracy. The third term on the right hand side of (13) can be calculated purely 

numerically. Thus we obtain the coefficients in (12): 

CgJ = -0.20665154, 

eg = -0.063212891, 

eg = -0.024202688, 

CgJ = -0.01062245, 

eg = -0.0050245436, 

eg = -0.002441955, 

eg = -0.0011555635. 

On the basis of these numbers one can obtain the following approximation: 

42)(^)^ 3 ^ 



-0.826606 + 0.179208 + 0.5387 - 0.350049 cj^ + 0.04253260^^ 
^ 1 + 0.559635^- 0. 196815 ' ^ ^ 

It should be noted that the purely gluonic cut in this diagram is zero according to the 
Landau- Yang theorem [TT] . 

(2) 

In the last step we take the imaginary part R^ , 6 = NA, A, S of the (10), (11), (14) 
(/i^ = m^): 

Rf =47rlmnf\uj), (15) 

1 + Wz-l 

uj = , . 

1-iy/z^ 
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The nonperturbative part of the polarization operator was calculated up to the oper- 
ators of dimension 8 in [12] . For the correlator of heavy quarks there are: 
the only operator of dimension 4, 

two operators of dimension 6, 

Ol = {g'r'^G^^Gl^G^J , Ol = (g'jX) , 
and seven operators of dimension 8: 

Ol = {{g'd^'^G^G'^^.f + l{g'Gl^G';,f) , Ol = {{gT'^G^G:,.)') , 

3 2 b b 2^2 2 4 2 b b 2 

O4 = {{g (T "G^^C^^) + -{g G^^^C^^) ) , 0^ = {{g f "G^^G"^^) ) , 

Ol = {{g' r'^GUljl) , Ol = {{g'r'^G^^Gl^G^^^.^,,) , Ol = {g'jX-J ■ 
Here j^" is the color current of the light quarks, gj^ = 0'^^,^.^^ = f ^ The 

q=u,d,s 

operator product expansion part Yl'^^^ of the polarization operator has the form ^] 
{y = <5^/(4m^), 2Fi{a,b,c, z) is hypergeometric function): 

+ 5^ g (f + g 4,. . (1. 2 + . 9/2. + (16) 

In this equation 

4o = -3/5, 4, = 0, ^2 = 0, ^3 = 24/7, = -116/7, 45 = 552/35, 



4o = 36/5, = 8/105, ^2 = 8/35, ^3 = 144/35, 4^ = 32/21, ^5 = -736/35, 

^4,1 = 








1 


2 


3 


4 


5 


6 


7 


1 





2/315 


-16/63 


92/21 


544/63 


-562/9 


160/3 





2 


9/35 


1/105 


32/63 


-118/9 


200/3 


-5639/45 


30448/315 


-180/7 


3 


-18/5 


-8/315 


64/63 


-752/21 


9920/63 


-10036/45 


320/3 





4 


-144/35 


-2/315 


-32/45 


-188/63 


7376/63 


-10586/45 


-26464/315 


1656/7 


5 


612/35 


-4/105 


-64/63 


184/7 


-4000/21 


20524/45 


-1984/7 


-432/7 


6 


72/35 





-128/315 


944/63 


-432/7 


-1544/45 


98944/315 


-1728/7 


7 


-324/35 


4/315 





40/9 


-1696/63 


292/3 


-14528/63 


1296/7 
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The as-correction of the first order to the gluon condensate was obtained in JH] • The 
corresponding term in the polarization operator has the following form: 



where P^{z) is given in equation (8) of [TB] . 



3 Moments and mass redefinition 

In our approach the moments contain the perturbative part up to terms, the oper- 
ator product expansion part with the operators of 4, 6 and 8 dimensions and perturbative 
correction of the first order to the lowest dimension operator: 

M„(g2) = j2 Mi'\Q')a\m') + M^^^ + 02M^'^'\Q')a{m') , (18) 

k=0 

The explicit formula for Mj^^ can be obtained analytically: 

= u ^Ji^l"":^^' 1 + ^' 5/2 + n, -y) , 

(4m^)" 41 [n + |) 

where y = Q'^/{Am?). The next perturbative terms in (18) we calculate numerically (see 
(3)): 

M^rO^) = ^ 7 R^'\r,m')dr 
" ^ (4m2)" 7 (r + ?/)"+i ■ 

Here i?W is given by (6), for the components of i?^^) see (7), (8), (9) and (15), (10), (11), (14). 
The operator product expansion moments can be easily obtained from (16): 

<"-2r(4;^(^ 

MO^'^ = ^ (<^+Vc4, ^,^ + \|^^'^/^^/^h Fi(l+n,3+^+n,ll/2+n,-y^^ . 

27(4m2)"+4V^n+i ^ ^'*(2 + i)!r(ll/2 + n) ' . / ' y;; 

As for a .(-correction to the operator O2, the corresponding moments can be obtained 
by numerical differentiation of (17): 



(4m2)2n!' dQ^' 8y^{l + y) ' 
8 



y = Q^/iAm"). 

Calculating the numerical values of the moments one can easily find that the a^- 
corrections to the moments are very large and in fact the series (18) is divergent. The 
traditional solution of this problem is the mass redefinition. Indeed, the pole mass in 
the above formulas has the meaning of the mass of free quark and is ill defined quantity 
in the charmonium sum rules. This problem was discussed in details in [3]. Following 

we consider MS scheme. The corresponding mass fh is taken at the scale /^^ = m^: 
fh = ffi{fn^). 

The pole mass m can be expressed in terms of fh as the perturbative series: 



2 — 2 

m = m 



k=l,2,. 

where Ki, K2, were obtained in [T^ : 

fsTi = 8/3 , K2 = 22.4162 , = 260.526 . 

(These numbers are given for ri; = 3.) 

Now we reexpand the series (18) in terms of the mass fh: 

MniQ') = jzM^n\Q^)o^\rn^) + M'''''' + 02M^^^^\Q^)a{m') . (19) 

A;=0 

Here 

M(i)(Q2) = M(i) - irinM(°) +K,{n + 1) Q^M^^l , 

(g^) = Mf - i^mM^i) + K,{n + 1) Q'm';^1 + n{^{n + 1) - K,)m^^^ + 

+ {n + 1)[K2 - Kl{n + 1))Q'm'^1 + ^{n + l){n + 2) Q^M^l , 
M02{i)(g2) ^ ^02{i) _ ^^(^ ^ 2)M„^2 + K,{n + 1) Q'M^^, . 

All moments in the right hand sides of these equations are computed with the MS mass 
fh. 

At the some other scale fx^ the function M^^ changes 

(Q2)^Mf(Q2) +M«(Q^)/5o In ^, a{m') ^ a{ix') (20) 



m 



to ensure the scale independence at order a^. 
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4 Restrictions on the gluon condensate value 



To analyse the obtained data we introduce the dimensionless ratio r„ of the moments 
(19): 



Here 5 stands for the continuum contribution. At large n 5 tends to zero: n —>■ oo, 6^0. 
The theoretical ratio of the moments depends on Q^, quark mass m, QCD coupling and 
condensates. First of all, let us fix Q^. At = the perturbative corrections as well as 
the higher terms of the operator product expansion series are very large. On the other 
hand, at large Q^, Q'^/{Afh?) > 3, the effective expansion parameter ajSolnlQ"^ /rh"^) in 
(20) becomes large. In papers the values Q'^/^im^) = 1, (5^/(4m2) = 2 were used. 

In the present paper we work at Q'^/{4fh'^) = 2. 

As for the QCD coupling constant, from the hadronic r-decay we know ^ 



as{ml) = 0.33 ±0.03, (22) 

rrir = 1.777 GeV is the mass of the r-lepton. To obtain as at any other scale we solve 
numerically the renormalization group equation. The choice of the scale is discussed in 
i: 

;x2 = Q2 + ^2 _ ^23) 

The value of the c-quark mass m is determined in 0,13] with high accuracy: 

rh = 1.275 ±0.015 GeV 

(this value is taken from [3 J. Now theoretical r„ depends on the condensates only. 

Using the value = 3.51051 ± 0.00012 GeV jS] we find the phenomenological esti- 
mation of Tn in (21): 

r„= ±(5 = 3.90 ±0.05 ±(5. (24) 

4m2 

The uncertainty appears mainly due to the error in the quark mass value. 
Now let us look on the table of the moments. 



n 




a2M(2)/M(o) 




aM02(i)/M02 


4 


0.097 


0.032 


0.083 


-0.83 


5 


0.070 


0.018 


0.20 


-0.51 


6 


0.041 


0.0013 


0.44 


-0.27 


7 


0.0087 


-0.017 


0.91 


-0.065 


8 


-0.026 


-0.037 


1.80 


0.12 


9 


-0.061 


-0.058 


3.47 


0.30 



In this table Q2/(4^2) ^ 2, (^G2)/(4m2) = 0.0002 and scale (23) are chosen. 
Operators Of, O4 are excluded. 

The values in all the columns in the certain line should be small enough (< 0.5 in 
magnitude), besides, the values in the third column should be smaller than in the second 
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one. All these requirements ensure the convergence of the series (19). Thus we see that 
there is very narrow range of n: n = 5 — 6, and one can construct the only ratio (21): 

Ms 
4m2M6 ■ 

For Q'^/[4:ffi^) = 1 there are no appropriate n at all. 

The values of the higher order operators, Of, O4, can be estimated with the help 
of factorization hypothesis or instanton model (for details see 0). In the instanton 
consideration 

1 2 

01 = ^0., O| = 

(ol,...,Ol) = (4,8,3,4,0,8,0)i^^O2. 

' r c 

Here instanton radius pc = 2.5GeV~^. The instanton concentration uq is connected with 
the gluon condensate: (^G^) = 32no. 

In the frame of the factorization hypothesis 

K--oO-(^4-i-i^-°-5-°)(o^)^- 

However, at n = 5 — 6 both models give negligible corrections. 

At first we put 5 = in (24). Then 3.85 < rg < 3.95. Now we find the interval of the 
gluon condensate values, where these restrictions hold: 

0.005 GeV^ < (— G^) < 0.006 GeV^ . 

vr 

The uncertainty in the coupling constant and variation of the scale (/x^ = or 
= + 2m^) give negligible correction to the value of gluon condensate. 
Now let us try to evaluate 6. Unlike to the vector channel, the mass of the second 
resonance with J^'-^ = l^"*" is unknown. If we suppose that the differences between 
two lowest resonances in vector and axial-vector channels are approximately equal (about 
0.6 GeV), the continuum threshold can be evaluated: sq ~ 4.1 GeV. Now one can compare 
the integrals like (3) in the intervals [4m^, 00) and [sq, 00) to evaluate the continuum 
contribution to the certain moment. In such way we find: 

5 < 0.2. 

It is especially important that 6 is positive, 6 > 0. 

It turns out that in the considered sum rules the increase of rg results in decrease of 
the value (— G^) (see fig.l). That is why introducing the positive 6 in (24) we obtain the 
smaller values of the gluon condensate. 

In some papers in the sum rules for heavy quarks the Coulomb-like corrections are 
summed up. This is legitimate way for the nonrelativistic problems only. However, in 
our case the quark velocities v ~ y^(1 -|- Q"^/ [Afh'^))/n ^0.7 are large enough, and the 
nonrelativistic corrections are not dominating. Detailed consideration of this question can 
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be found in [Hj. 

Our final result is: 



i—G^) = (0.005 + 0.001 - 0.004) GeV^ . (25) 

71 

Note again that the continuum (i.e. higher states) does not affect the upper limit in (25), 
which can be considered as quite reliable. As for the lower limit, it depends on 6, for 
which one has the estimation only. Therefore, even zero value of the condensate can not 
be certainly excluded. 

Our result (25) is in agreement with the values of the gluon condensate, obtained 
in P] ii^G^) = 0.009 ± 0.007 GeV^) and in ^5J ((^G^) < 0.008 GeV^). However, it is 
significantly smaller than 0.012 GeV^, originally obtained in PP. 

The author is thankful to B.L.Ioffe for the formulation of the problem and fruit- 
ful discussions and to K.N. Zyablyuk for helpful discussions. 

The work is supported in part by grants INTAS 2000 Project 587 and RFFI 03-02- 
16209. 
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Figure 1: The ratio rg as the function of the gluon condensate. Horizontal hnes denote the 
hmits in (24). 

References 

[1] M. Shifman, A. Vainshtein and V. Zakharov, Nucl.Phys. B147 (1979) 385, 448. 
[2] S.Nikolaev and A. Radyushkin, JETP Lett. 37 (1982) 526. 
[3] B. loffe and K. Zyablyuk, Eur.Phys. J. C27 (2003) 229. 



12 



[4] B. loffe, Phys.Atom.Nucl. 66 (2003) 30; Yad.Fiz. 66 (2003) 32. 

[5] K.Zyablyuk, JHEP 0301 (2003) 081. 

[6] K. Hagiwara et al. (Particle Data Group), Phys. Rev. D66 (2002) 010001. 

[7] A.H. Hoang, T.Teubner, Nucl.Phys. B519 (1998) 285. 

[8] K.Chetyrkin, J. Kuhn, M. Steinhauser, Nucl.Phys. B505 (1997) 40. 

[9] K. Chetyrkin, R. Harlander, M. Steinhauser, Phys.Rev. D58 (1998) 014012. 

[10] B. Kniehl, J. Kuhn, Nucl.Phys. B329 (1990) 547. 

[11] L. Landau, Dokl.Akad. Nauk USSR 60 (1948) 207; C.Yang, Phys.Rev. 77 (1950) 242. 

[12] S.Nikolaev and A. Radyushkin, Sov.J.Nucl.Phys. 39 (1984) 91; 
Yad.Fiz. 39 (1984) 147. 

[13] D. Broadhurst, P. Baikov, V. Ilyin, J. Fleischer, O.Tarasov, V. Smirnov, Phys. Lett. 
B329 (1994) 103. 

[14] N.Gray, D. Broadhurst, W. Grafe, and K. Schilcher, Z.Phys. C48 (1990) 573; 
K. Melnikov and T. vanRitbergen, Phys.Lett. B482 (2000) 99; K.Chetyrkin and 
M. Steinhauser, Nucl.Phys. B573 (2000) 617. 

[15] B. Geshkenbein, B. loffe and K.Zyablyuk, Phys.Rev. D64 (2001) 093009. 



13 



